Multi-scale modeling and variability in cardiac cellular electrophysiology

Cardiac rhythm is the result of interactions between millions of heart muscle cells. Genetic, molecular and functional factors are involved on the cellular, tissue and organ level. Multi-scale computer models are used to study this complex interplay of processes on various levels. This dissertation introduces a software application that makes these models more versatile and easier to use. It also presents results on the effects of genetic mutations and natural variability on ion channel function in heart muscle cells. These results improve our understanding of the mechanisms of the heart and can serve as a catalyst for studies on cardiac arrhythmias

Physiology-based regularization of the electrocardiographic inverse problem

The inverse problem of electrocardiography aims at noninvasively reconstructing electrical activity of the heart from recorded body-surface electrocardiograms. A crucial step is regularization, which deals with ill-posedness of the problem by imposing constraints on the possible solutions. We developed a regularization method that includes electrophysiological input. Body-surface potentials are recorded and a computed tomography scan is performed to obtain the torso-heart geometry. Propagating waveforms originating from several positions at the heart are simulated and used to generate a set of basis vectors representing spatial distributions of potentials on the heart surface. The real heart-surface potentials are then reconstructed from the recorded body-surface potentials by finding a sparse representation in terms of this basis. This method, which we named 'physiology-based regularization' (PBR), was compared to traditional Tikhonov regularization and validated using in vivo recordings in dogs. PBR recovered details of heart-surface electrograms that were lost with traditional regularization, attained higher correlation coefficients and led to improved estimation of recovery times. The best results were obtained by including approximate knowledge about the beat origin in the PBR basis

Rapid Characterization of hERG Channel Kinetics I: Using an Automated High-Throughput System

Predicting how pharmaceuticals may affect heart rhythm is a crucial step in drug-development, and requires a deep understanding of a compound’s action on ion channels. In vitro hERG-channel current recordings are an important step in evaluating the pro-arrhythmic potential of small molecules, and are now routinely performed using automated high-throughput patch clamp platforms. These machines can execute traditional voltage clamp protocols aimed at specific gating processes, but the array of protocols needed to fully characterise a current is typically too long to be applied in a single cell. Shorter high-information protocols have recently been introduced which have this capability, but they are not typically compatible with high-throughput platforms. We present a new 15 second protocol to characterise hERG (Kv11.1) kinetics, suitable for both manual and high-throughput systems. We demonstrate its use on the Nanion SyncroPatch 384PE, a 384 well automated patch clamp platform, by applying it to CHO cells stably expressing hERG1a. From these recordings we construct 124 cell-specific variants/parameterisations of a hERG model at 25C. A further 8 independent protocols are run in each cell, and are used to validate the model predictions. We then combine the experimental recordings using a hierarchical Bayesian model, which we use to quantify the uncertainty in the model parameters, and their variability from cell to cell, which we use to suggest reasons for the variability. This study demonstrates a robust method to measure and quantify uncertainty, and shows that it is possible and practical to use high-throughput systems to capture full hERG channel kinetics quantitatively and rapidly

Reproducible model development in the Cardiac Electrophysiology Web Lab

The modelling of the electrophysiology of cardiac cells is one of the most mature areas of systems biology. This extended concentration of research effort brings with it new challenges, foremost among which is that of choosing which of these models is most suitable for addressing a particular scientific question. In a previous paper, we presented our initial work in developing an online resource for the characterisation and comparison of electrophysiological cell models in a wide range of experimental scenarios. In that work, we described how we had developed a novel protocol language that allowed us to separate the details of the mathematical model (the majority of cardiac cell models take the form of ordinary differential equations) from the experimental protocol being simulated. We developed a fully-open online repository (which we termed the Cardiac Electrophysiology Web Lab) which allows users to store and compare the results of applying the same experimental protocol to competing models. In the current paper we describe the most recent and planned extensions of this work, focused on supporting the process of model building from experimental data. We outline the necessary work to develop a machine-readable language to describe the process of inferring parameters from wet lab datasets, and illustrate our approach through a detailed example of fitting a model of the hERG channel using experimental data. We conclude by discussing the future challenges in making further progress in this domain towards our goal of facilitating a fully reproducible approach to the development of cardiac cell models

Models of the cardiac L-type calcium current: A quantitative review

The L-type calcium current (ICaL) plays a critical role in cardiac electrophysiology, and models of ICaL are vital tools to predict arrhythmogenicity of drugs and mutations. Five decades of measuring and modeling ICaL have resulted in several competing theories (encoded in mathematical equations). However, the introduction of new models has not typically been accompanied by a data-driven critical comparison with previous work, so that it is unclear which model is best suited for any particular application. In this review, we describe and compare 73 published mammalian ICaL models and use simulated experiments to show that there is a large variability in their predictions, which is not substantially diminished when grouping by species or other categories. We provide model code for 60 models, list major data sources, and discuss experimental and modeling work that will be required to reduce this huge list of competing theories and ultimately develop a community consensus model of ICaL. // This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Molecular and Cellular Physiolog

Autocorrelated measurement processes and inference for ordinary differential equation models of biological systems

Ordinary differential equation models are used to describe dynamic processes across biology. To perform likelihood-based parameter inference on these models, it is necessary to specify a statistical process representing the contribution of factors not explicitly included in the mathematical model. For this, independent Gaussian noise is commonly chosen, with its use so widespread that researchers typically provide no explicit justification for this choice. This noise model assumes `random' latent factors affect the system in ephemeral fashion resulting in unsystematic deviation of observables from their modelled counterparts. However, like the deterministically modelled parts of a system, these latent factors can have persistent effects on observables. Here, we use experimental data from dynamical systems drawn from cardiac physiology and electrochemistry to demonstrate that highly persistent differences between observations and modelled quantities can occur. Considering the case when persistent noise arises due only to measurement imperfections, we use the Fisher information matrix to quantify how uncertainty in parameter estimates is artificially reduced when erroneously assuming independent noise. We present a workflow to diagnose persistent noise from model fits and describe how to remodel accounting for correlated errors

Accounting for variability in ion current recordings using a mathematical model of artefacts in voltage-clamp experiments

Mathematical models of ion channels, which constitute indispensable components of action potential models, are commonly constructed by fitting to whole-cell patch-clamp data. In a previous study, we fitted cell-specific models to hERG1a (Kv11.1) recordings simultaneously measured using an automated high-throughput system, and studied cell-cell variability by inspecting the resulting model parameters. However, the origin of the observed variability was not identified. Here, we study the source of variability by constructing a model that describes not just ion current dynamics, but the entire voltage-clamp experiment. The experimental artefact components of the model include: series resistance, membrane and pipette capacitance, voltage offsets, imperfect compensations made by the amplifier for these phenomena, and leak current. In this model, variability in the observations can be explained by either cell properties, measurement artefacts, or both. Remarkably, by assuming that variability arises exclusively from measurement artefacts, it is possible to explain a larger amount of the observed variability than when assuming cell-specific ion current kinetics. This assumption also leads to a smaller number of model parameters. This result suggests that most of the observed variability in patch-clamp data measured under the same conditions is caused by experimental artefacts, and hence can be compensated for in post-processing by using our model for the patch-clamp experiment. This study has implications for the question of the extent to which cell-cell variability in ion channel kinetics exists, and opens up routes for better correction of artefacts in patch-clamp data. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'

Model-driven optimal experimental design for calibrating cardiac electrophysiology models

Background and Objective: Models of the cardiomyocyte action potential have contributed immensely to the understanding of heart function, pathophysiology, and the origin of heart rhythm disturbances. However, action potential models are highly nonlinear, making them difficult to parameterise and limiting to describing ‘average cell’ dynamics, when cell-specific models would be ideal to uncover inter-cell variability but are too experimentally challenging to be achieved. Here, we focus on automatically designing experimental protocols that allow us to better identify cell-specific maximum conductance values for each major current type.Methods and Results: We developed an approach that applies optimal experimental designs to patch-clamp experiments, including both voltage-clamp and current-clamp experiments. We assessed the models calibrated to these new optimal designs by comparing them to the models calibrated to some of the commonly used designs in the literature. We showed that optimal designs are not only overall shorter in duration but also able to perform better than many of the existing experiment designs in terms of identifying model parameters and hence model predictive power.Conclusions: For cardiac cellular electrophysiology, this approach will allow researchers to define their hypothesis of the dynamics of the system and automatically design experimental protocols that will result in theoretically optimal designs

Immediate and Delayed Response of Simulated Human Atrial Myocytes to Clinically-Relevant Hypokalemia

Although plasma electrolyte levels are quickly and precisely regulated in the mammalian cardiovascular system, even small transient changes in K+, Na+, Ca2+, and/or Mg2+ can significantly alter physiological responses in the heart, blood vessels, and intrinsic (intracardiac) autonomic nervous system. We have used mathematical models of the human atrial action potential (AP) to explore the electrophysiological mechanisms that underlie changes in resting potential (Vr) and the AP following decreases in plasma K+, [K+]o, that were selected to mimic clinical hypokalemia. Such changes may be associated with arrhythmias and are commonly encountered in patients (i) in therapy for hypertension and heart failure; (ii) undergoing renal dialysis; (iii) with any disease with acid-base imbalance; or (iv) post-operatively. Our study emphasizes clinically-relevant hypokalemic conditions, corresponding to [K+]o reductions of approximately 1.5 mM from the normal value of 4 to 4.5 mM. We show how the resulting electrophysiological responses in human atrial myocytes progress within two distinct time frames: (i) Immediately after [K+]o is reduced, the K+-sensing mechanism of the background inward rectifier current (IK1) responds. Specifically, its highly non-linear current-voltage relationship changes significantly as judged by the voltage dependence of its region of outward current. This rapidly alters, and sometimes even depolarizes, Vr and can also markedly prolong the final repolarization phase of the AP, thus modulating excitability and refractoriness. (ii) A second much slower electrophysiological response (developing 5–10 minutes after [K+]o is reduced) results from alterations in the intracellular electrolyte balance. A progressive shift in intracellular [Na+]i causes a change in the outward electrogenic current generated by the Na+/K+ pump, thereby modifying Vr and AP repolarization and changing the human atrial electrophysiological substrate. In this study, these two effects were investigated quantitatively, using seven published models of the human atrial AP. This highlighted the important role of IK1 rectification when analyzing both the mechanisms by which [K+]o regulates Vr and how the AP waveform may contribute to “trigger” mechanisms within the proarrhythmic substrate. Our simulations complement and extend previous studies aimed at understanding key factors by which decreases in [K+]o can produce effects that are known to promote atrial arrhythmias in human hearts